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Abstract 

We apply the soft-collinear effective theory (SCET) to deep inelastic scattering near the endpoint 
region. The forward scattering amplitude, and the structure functions are shown to factorize as 
a convolution of the Wilson coefficients, the jet functions, the parton distribution functions. The 
behavior of the parton distribution functions near the endpoint region is considered. It turns 
out that it evolves with the Altarelli-Parisi kernel even in the endpoint region, and the parton 
distribution function can be factorized further into a collinear part and the soft Wilson line. The 
factorized form for the structure functions is obtained by the two-step matching, and the radiative 
corrections or the evolution for each factorized part can be computed in perturbation theory. We 
present the radiative corrections of each factorized part to leading order in a^, including the zero-bin 
subtraction for the collinear part. 
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I. INTRODUCTION 

The soft-collinear effective theory (SCET) is a useful theoretical tool to treat 

physical processes with energetic light particles in a systematic way. For an energetic particle 
moving in the n'^ direction, the momentum can be decomposed into 

= ~ 0{ff) + 0(A) + O(AVQ), (1) 

where Q is a large scale, and ra^, W- are lightlike vectors satisfying ra^ = = 0, n ■ 7T = 2. 
Each component has a distinct scale in powers of A which is a typical hadronic scale, and 
SCET describes the interactions of the collinear particles and the ultrasoft (usoft) particles 
with momentum ~ (A, A, A). Since there are three distinct scales for the momentum 
of a collinear particle, SCET employs a two-step matching process by integrating out large 



energy scales successively js]. In the first stage the degrees of freedom of order Q from 
the full theory are integrated out to produce SCETj. In SCETi, collinear particles are 
allowed to interact with usoft particles and the typical virtuality of the collinear particles 
is ~ QK. In the second stage the degrees of freedom with ~ QK are integrated out, 
and the remaining effective theory in which all the particles have ~ A^ is called SCETn- 
Here the collinear particles are decoupled from the soft particles. The Wilson coefficients 
of operators and the renormalization behavior of them can be computed perturbatively by 
matching the effective theories at each boundary. 

SCET has been successfully applied to various B meson decays 



It is especially convenient to study the factorization properties of B decays including 



spectator interactions since SCET is formulated such that soft and collinear particles are 



decoupled. On the other hand, SCE' 
include energetic light particles 



can b e applie d to other high-energy processes which 
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DIS) near the endpoint region using SCET Il4,ll7| and using an effective theory 



)een applied to deep inelastic 



In this paper we analyze the endpoint region in DIS more carefully using the two-step 
matching to show the explicit factorization of the structure functions in terms of the hard 
part, the jet function, the soft gluon emissions, and the collinear matrix elements. We also 
discuss and compare delicate physical meanings and implications of the parton distribution 
functions in the endpoint region, defined both in the full theory and in SCET. In addition 
to showing the factorization, we take one step further to consider another aspect of DIS, 



namely the behavior of the longitudinal structure function near the endpoint region. The 

longitudinal structure function vanishes at leading order in as due to the fact that the parton 

(quark) in the proton has spin 1/2. However this is broken at order and the longitudinal 

structure function is further suppressed by A/Q, which we explicitly present here. 

In Section [n] we explain the kinematics of DIS. We choose the Breit frame and present 

how the momenta scale in powers of A, which is useful in constructing and matching effective 

theories. The forward scattering amplitude, and the structure functions are defined in SCET, 

compared with those in the full theory. In Section IIIII the method to compute the forward 

scattering amplitudes in DIS using SCET is described. The leading and the subleading 

currents are introduced and the prescription for the usoft factorization is explained. In 

Section HVl we compute the structure function Fi{x,Q), and show that it factorizes. In 

Section |Vl we consider the parton distribution near the endpoint region, and express the 

forward scattering amplitude in terms of the parton distribution function. In Section IVIl 

we present the moments of the structure functions as a product of the moments for each 

factorized term. In Section IVlIl we compute the radiative corrections of each factorized term 

to order in a^, and express the moments of the structure functions to leading logarithmic 

accuracy. In Section IVIIIl we compute the longitudinal structure function Fl{x, Q) in SCET 

and show that it also factorizes. In the final section, we give a conclusion. In Appendix A, 

□ 

the zero-bin subtraction method [19| in SCETi before the usoft factorization of the collinear 
fields is discussed. In Appendix B, the procedure for taking the imaginary part in SCETi 
and SCETn is explained. In Appendix C, the anomalous dimension of the operator J^^''^ is 
computed to order a^. 

II. KINEMATICS 

Let us consider the electroproduction in DIS ep — >■ eX near the endpoint region. The 
hadronic process consists of 7*p — ^ X, and we choose the Breit frame in which the incoming 
proton is in the direction, and the outgoing hadrons are mainly in the direction. The 
momentum transfer from the leptonic system is given by 

q^ = {ji- g, g^:, n ■ q) = (Q, 0, -Q) = |k - n^, (2) 
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where — —Q"^ is the large scale. The Bjorken variable x is defined as 

Q ... 

where is the proton momentum in the direction. The momenta of the proton P^ and 

the final-state particles px = P + q are given by 

P^ = (rr-P,Pj:,n-P)~(^,Pi:,|), 

_ 1 — X 

Px = [ri-px:Px±:n-px 



X 

32 _ , A 2 ^2 _ /n2 



with P ~ A , = g (1 — x)/x where A is a typical hadronic scale of order 1 GeV. Near 
the endpoint where x approaches 1 (1 — x ~ A/g)^, the invariant mass squared of the 
final-state particles becomes ~ g^(l — x) ~ gA. Then the final-state particles can be 
regarded as coUinear particles in SCETi, which are integrated out to obtain SCETn through 
the two-step matching procedure. 

At the parton level, let p^ be the momentum of the incoming parton inside the proton, 
and let y be the longitudinal momentum fraction [n-p = yn- P). Then the partonic Bjorken 
variable w is given as 

Q'^ n-q X 

w — ~ = — . (oj 

2p ■ q n-p y 

The momentum of the outgoing parton as p'^ can be written as 

p'M ^ + ^ . p'^^^. ^ (g^ p>^^ {l-w)n-p). (6) 

And the endpoint region corresponds to 1 — w ^ K/Q such that p'^ ~ gA. 
The spin-averaged cross section for DIS can be written as 

where k and k' are the incoming and outgoing lepton momenta with q — k' — k, L^^ is the 
lepton tensor, and s = (p + k)'^. The hadronic tensor W^j, is related to the imaginary part of 
the forward scattering amplitude T'^'^ . The forward scattering amplitude is the spin- averaged 
matrix element of the time-ordered product of the electromagnetic currents, written as 

T^.{x,Q) = (P|f^.|P)spin av., %.{x,Q) = I J dSe^^"T[jt (^) J,(0)] , (8) 



^ In fact, 1 — X does not have to be of order A/Q. Instead, we can introduce a small parameter 6 = 1 — x 
and Px ^ Q5 » A^. But for simplicity, we consider the case with 1 — a; ~ A/Q. 
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where is the electromagnetic current. The relation between the hadronic tensor W^i, and 
the forward scattering amplitude T^u is given by 

W^,ix,Q) = -lmT^,ix,Q). (9) 

TT 

In electroproduction, considering all the possible Lorentz structure, T^^, can be generally 
written as 

T^u{x, Q) = -gf^JTi + in^n^ + n^n^)T2 + {n^n^ - n^n^)T3 + n^n^Ti + n^n^T^, (10) 

where gj^^^ = g^y — {n^riu + ?T^?2,^)/2. Due to the current conservation {q^T'^^ = 0), and the 
parity conservation, we have T4 = T5 = T2, and T3 = 0. Therefore the forward scattering 
amplitude has two independent quantities, and is given by 

Ti„y{x, Q) = -g^^Ti{x, Q) + {n^ + + n^)T2{x, Q). (11) 

This can be cast into different forms using the fact that and any terms proportional to 
Q/i = Q{n^ — n^)/2 can be discarded since they vanish when they are contracted with the 
lepton tensor. We can write + rl^ = (rT^ — n^) + 2n^ = (ra^ — n^) + 2n^ and drop the 
terms proportional to ri^ — ?T^. Then Eq. ( TTTl) can be equivalently written as 

T^,(x, Q) = -gfMx, Q) + 4n^n,T2(x, Q) = -g^Ji{x, Q) + m^n^T^^x, Q). (12) 

The structure functions are defined from the hadronic tensor Wf^i, as 

W^,{x, Q) = -g,.F,{x, Q) + ^H^, Q), (13) 

where the terms proportional to or q^, are dropped. Using = n ■ 2P ■ q = 

n ■ Pn ■ q = Q'^/x, we can write Eq. (fT3l) as 

77 77 

W^^ix.Q) = -g^,F^{x,Q) + ^F2{x,Q) 

2x 



1 Tl Tl 

tT tT /I \ 

-g;^^F^{x,Q) + ^{-F2{x,Q) - 2F,{x,Q)^ 
= -giMx,Q) + ^FL{x,Q), (14) 

where we extract — 77^ and discard it to obtain the third relation, and the longitudinal 
structure function Fl{x,Q) is defined as 

Fl{x,Q) = -F2{x,Q)-2F,{x,Q). (15) 



The Lorentz structure n^ny in the final expression of Eq. (fT^ can be replaced by n^riy. 
Comparing Eqs. ffTTj) and (|T^ . we obtain the relations 

Q) = -lmTi(x, Q), Fl{x, Q) = -lmT2(x, Q). (16) 

TT TT 

As we will show explicitly, Ti(a;, Q) receives the contribution at leading order, and T2{x, Q) 
is suppressed by A/Q and compared to Ti{x,Q). Therefore the Callan-Gross relation 

= holds to leading order, but is violated at subleading order. Here we also present 
computed using SCET. In fact, the equivalence between Eq. (ITT]) and Eq. (fT2!) turns out 
to imply nontrivial relations because the subleading contributions proportional to n^n'^ and 
Y[^^Y[^J■ come from different subleading current operators in SCET. The statement that all the 
expressions are equivalent means that the longitudinal structure functions can be obtained 
using any subleading current operators and it holds to all orders in a^. The nontrivial 
relation will be verified in this paper at order a^. 



III. OPERATORS IN SCET NEAR THE ENDPOINT REGION 



In computing the forward scattering amplitude, we first express the electromagnetic cur- 
rent in terms of the effective fields in SCETj. The electromagnetic current operator at 
leading order in SCET is given by 



(17) 



where ^„ {^n) is the n (n) coUinear fermion field in SCET. Here Wn and Wn are the coUinear 
Wilson lines 



Wjx) 



[ exp(-5f^^?2 ■ A„(x))] , Wn{x) = [ exp(-5(;3^^?2 ■ An{x) 

perms ^ ' ' perms 



:i8) 



It. / perms ' 

Here Al^ i^n) is the coUinear gluon in the {n'^) direction and the summation over the 
label momenta is suppressed. The Wilson coefficient C{Q,^) is actually an operator and 
Eq. f|T7|) is written as 



= J dooduCiuj, Lj)'^^WJ{uj - n ■ V^h^diuj - n ■ V)Wl + h-c, (19) 

where n-V in ■ V'^) is the operator extracting the label momentum in the n (n) direction. 
The operator form in Eq. (|T9|l is useful in deriving the Feynman rules to compute radiative 



corrections. The hard coefficient C{Q,ji) can be obtained from matching the full theory 
onto SCETi, it is given to order as as 

C(g.,)^l + ^(-.„^f + 3.„f-84). (20, 

The hard coefficient C{Q,fi) satisfies the renormalization group equation 

A^^^=7HMC(Q,ri. 7«M = ^^^(41u^ + 3). (21) 



We can obtain subleading current operators at order yA/Q, which contain either iD:^ 
or iD^. There are two independent operators involving iD:^, one of which arises from the 
subleading correction to the fermion field C,n 

11, q - LWn{i + ^wl^w^^ + ■ ■ ^iMU- (22) 

The second term in Eq. ( 122|) yields the subleading current at tree level 

^i"^^ = L^^W^^^.Wl^n = -n.lMWn^^WlU (23) 

The second type arises from integrating out the off-shell modes when the coUinear quark 
emits a coUinear gluon A^, and it is given at tree level as 

Jll'^ = lnl,^nWn\-^Wlin = -n^l^Up^W^^Wlin. (24) 

This can be derived by computing the Feynman diagram for the process and by integrating 
out the intermediate state with virtuahty ~ Q^. And the result can be made gauge 
invariant by inserting the appropriate coUinear Wilson lines. A novel method to derive the 
operator is the auxiliary field method 



There are other subleading current operators involving ^-D^, which can be obtained by 
expanding to subleading order and by considering the process in which ^„ emits Af^. 
However these subleading operators do not contribute to the jet function which is obtained 
by integrating out the degrees of freedom of order ~ QK in going down to SCETn because 
these subleading operators describe the interaction of the n-coUinear particles in the proton. 
These operators contribute to the subleading corrections for the parton distribution functions 
which are given by the matrix elements of the coUinear operators in the direction, and 
we will not consider them here. 
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Before going down to SCETn, it is convenient to factor out the usoft interactions by 
redefining the colhnear fields, for example, as 



c p /t^Vf P » Y-P- Y- A>^Y^ 



(25) 



for the colhnear fields moving from — oo to x. Once the usoft interactions are factored out, 
collinear particles do not interact with usoft particles any more. The prescription of the 
usoft Wilson lines depends on the propagation of the collinear particles or antiparticles to 
which the soft gluons are attached, and it is described in detail in Ref. [l^. The possible 
usoft Wilson lines are given by 

1 



J2 

5erm 



perm 



1 



n -TZ — ie 



'-gn ■ Aus) 



[-gn ■ Ay,,) 



Yn{x) = Pexp 
Yn{x) = Pexp 



^9 



(26) 



where TZ is the momentum operator for the usoft fields and the path ordering P means that 
the fields are ordered such that the gauge fields closer (farther) to the point x are moved to 
the left, while P denotes the anti-path ordering. As explained in Ref. 16|, Yn iY^) is the 
usoft Wilson line attached to the collinear particle (antiparticle) from — cxd, while Yn iY^) 
is the usoft line attached to the collinear antiparticle (particle) moving to oo. This delicate 
procedure of choosing the appropriate usoft Wilson lines is related to the ie prescription, 
which specifies the location of the poles. Physically, this is related to choosing the sign of 
n ■ p since the denominator in Eq. ( l26l) is actually n -TZ + isgn.{n ■ p)e, and the sign oin-p 
depends on whether the collinear field is a particle or an antiparticle. 

Now that the current operators at leading and subleading order in SCET are known, we 
can compute the forward scattering amplitude T^^, or the hadronic tensor W^^, and factorize 
the usoft interactions using the appropriate prescription for the usoft Wilson lines. Then 
we integrate out the degrees of freedom of order p'^ ~ QA to obtain the result in SCETn. 
In SCETn the soft interactions are decoupled from the collinear particles with p"^ ~ A^, and 
the decoupled soft particles contribute to the soft Wilson lines which are responsible for the 
emission of soft gluons. 
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(b) 



— oo 

FIG. 1: (a) Feynman diagram for the forward scattering amplitude in DIS in SCETi and (b) the 
prescription of the (u)soft Wilson lines. 

IV. FACTORIZATION OF Fi(a;,Q) 

In SCETi after the usoft factorization, the time-ordered product Tj^^^ at leading order is 
written as 

= zC'iQ) J d'ze'^'^+^-^'>^T[^^WnYh,Y^WlU^)^^^^ (27) 

where and p''^ are the label momenta. Note that n-q = n-p' which ensures the conservation 
of the label momenta in the direction, while there is a slight mismatch in the direction 
near the endpoint w ~ 1 such that n- q + n-p = {l — w)n-p, which survives in the exponent. 
The Feynman diagram of the forward scattering amplitude for T^°) is sketched in Fig. [1] (a). 

The prescription for the usoft Wilson lines in DIS is described in Fig. [T] (b). It is deter- 
mined by the external states, which consist of an incoming particle C,n from — oo to 0, and 
an outgoing particle from z to oo. The intermediate states can move either from to 
— oo and then from — oo to z, or from to oo and then from oo to z. In both cases, the usoft 
Wilson line survives between and z, and the remaining part cancels. Either choice of the 
intermediate states is appropriate for describing DIS and here we choose the usoft Wilson 
lines for each current as 

InWnl.Wlir. : In ^ IM K ^ ^nA^M. in - Yni. ^ Y^A^^Y^ (28) 

where the intermediate state is going from to — oo, then moving from — oo to z. This 
prescription is used in Eq. (j27ll . 
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Since there are no collinear particles in the direction in the final state, we obtain the 
jet function defined by 



-ik-z 



Jp{k), 



(29) 



2 J (27r)4 

where P is the label momentum and Jp{k) depends only on n-k. We can simplify T^°^ using 
the fact that 

d^k 



d^z 



-ik±-z±—in-kn-z/2 



[ d'^z — / dn ■ k6'^(z±)6( ) = — [ dn ■ kdn ■ 

J Att J ^ ^ V 2 ^ 47r 7 



z, (30) 



and plugging the jet function into Eq. (1271) . we have 

f (?) = -C\Q) [ duo 
X T 



dn ■ Zdn ■ k^i^^/2+n-q-n-k]n-z/2j^^^ . 



4lT 

dn ■ zdn ■ k 



-C'iQ) / duj 



X T 



^^^i[u^/2+n-q-n-k-n]n-z/2 j^^^ . 



— C^iQ) J duj J dn ■ k J drjdi^— + n ■ q — n ■ k — 7]^Jp{n ■ k) 



X -(0|tr 



~SlSJ{r^ + n ■ id)SiSn] |0) i^WnS{uj - V+)^*-lMU^) 

ifi 



= ~C\Q) J duj dr^Jp(^ +n-q- 7^)S{v)UWnS{uJ - V+)^,^^M^n, (31) 

where the usoft Wilson line F„ (Yn) in SCETj is replaced by the soft Wilson line S'„ (Sn) in 
SCETn- Here the operator P4. = n ■ V + n ■ V'^ is the sum of the label momenta. Since the 
soft interaction is decoupled from the collinear sector, the soft Wilson lines are pulled out, 
and are described by the vacuum expectation of the soft Wilson line S{ri), which is given by 

S{v) = j^{0\tT[slSJ{v + n . zd)SiSn] |0). (32) 

The delta function in Eq. (!3T!) states that the momentum conservation in the n direction 
includes the soft momentum from soft gluons. Eq. fl3Tl) is the factorized form for the leading 
forward scattering amplitude. It consists of the hard part C^{Q), obtained in matching the 
current between the full theory and SCETj, the jet function Jp{n ■ k), obtained in matching 
between SCETj and SCETn, and the remaining collinear and soft operators in SCETn, 
whose matrix elements are given by nonperturbative parameters. The radiative corrections 
or the renormalization group evolution of each term can be computed in perturbation theory. 
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Note that the final operators in SCETn show a pecuhar structure. In inclusive B decays, 
the final operator after the two-step matching is a heavy quark bilinear operator with the soft 
Wilson lines. The matrix element of this operator is parameterized by the shape function of 
the B meson [20]. This is because the final operator is made of soft particles. But in DIS, the 
final operators are made of the collinear operators and the soft Wilson line. The soft Wilson 
line is responsible for the soft gluon emission and the parton distribution function near the 
endpoint is affected by this when a collinear particle participates in the hard scattering. It is 
a general feature for physical processes with collinear external particles that the soft Wilson 
line does not completely cancel near the endpoint, and it describes the soft gluon emission 
in the process. 

V. PARTON DISTRIBUTION FUNCTION 

We can extract Ti^\ proportional to —g^^, in T^^^ from Eq. fl3Tl) . and it is given by 



271 

X {P\^^Wn6{uj-V+)'^WlumP)s:>in..., (33) 

where z = n-z/2. We want to express Eq. ( I33l) in terms of the parton distribution functions. 
Here we consider only the flavor nonsinglet contribution. The standard coordinate space 



deflnitions 211] for the proton parton distribution functions fp{y) for quarks of flavor q 
moving in the n direction in full QCD are given as 

where Y{y, 0) is the path-ordered Wilson line and \P{P)) is the proton state with momentum 
P^. Here y is deflned as the longitudinal momentum fraction of the parton before the hard 
scattering and that of the proton, n ■ p = yn ■ P. 

The deflnition of the parton distribution function in Eq. flM|l is appropriate away from 
the endpoint region. But near the endpoint region, we have to extend the deflnition of the 
parton distribution to include the effect of the soft gluon emission, satisfying the requirement 
that it approach Eq. (1341) away from the endpoint region. At flrst sight, the soft momentum 
does not affect the parton distribution function since it describes the large energy component 
of the parton. In order to see why this is not so, let us consider a parton near the endpoint 
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region undergoing a hard collision, as depicted in Fig. [2l First a parton with the longitudinal 
momentum fraction y comes out of the proton. It emits soft gluons with total momentum 
K ~ Aqcd before it undergoes a hard collision with a photon. The parton distribution 
function fp{n-p/n ■ P) describes the probability of a parton entering the hard collision with 
the longitudinal momentum fraction n ■ p = yn ■ P — k. When yn ■ P ^ Aqcd including 
the endpoint region, the inclusion of k, seems to give a negligible effect. When we take a 
time-ordered product of this current as in Fig. [1] (a), all the soft gluons are attached to the 
n-collinear outgoing fermion due to the property of the soft interactions, which means that 
all the soft gluons are real gluons when we take the discontinuity. Away from the endpoint 
region, n ■ p' and n ■ p' of the n-coUinear quark are of order Q, and are not affected by the 
interaction of the soft gluons, that is, n-p' does not change to leading order in A. Therefore 
the interaction of the soft gluons can be neglected and we can safely put n = 0. 

Near the endpoint region, however, n ■ p' is of order A and the interaction with the soft 
gluons can significantly affect n -p'. If a physical quantity depends on the term proportional 
to 1/n-p', like the jet function, we have to keep the momentum k of the soft gluons carefully. 
From the above argument, the parton distribution function in Eq. ( IMl) can be extended in 
the endpoint region as 

fUy,^) = J ^e-(^-^-'')(P(P)|g(^)r(.,0)^g(0)|P(P))3pinav., (35) 

where the naive yn ■ P is replaced hj yn ■ P — k, that is, the parton distribution function 
is a function of the large longitudinal momentum fraction y, and the momentum of the soft 
gluons K. The Wilson line Y{z, 0) in Eq. is the Wilson line with the gauge field in 




FIG. 2: An energetic parton comes out of the proton with the momentum yn ■ P. It emits soft 
gluons with momentum k, and the momentum of the hard parton before the hard scattering 
becomes yn ■ P — k. 
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full QCD. However, by looking at the kinematics near the endpoint, the Wilson line can be 
decomposed into the coUinear and the soft Wilson lines [22 1. This procedure is similar to 
the approach in SCET, and SCET makes this procedure manifest. 
The parton distribution function in SCETn can be written as 

2tc J 



xSiv){Pn\[UWn]l6iu-P, 



(36) 



where the spin average is implied. The additional exponential factor e*'^^/^ comes from the 
label momentum of C,n- The usoft Wilson lines are prescribed according to Eq. (!28l) . and 
the proton state \P) is replaced by \Pn), in which the valence quarks are collinear in the n 
direction. 

Let us define a new parameter g'p{u), given by the spin-averaged matrix element of the 
collinear operators, as 



{Pn\UWn5{uj - V+)^Wlin\Pn) = n-P I du6{u - 2un ■ P)g%{u) 



(37) 



where the contribution from the antiquark is discarded for simplicity. Physically gp{u) cor- 
responds to the probability for the proton to emit a parton with the longitudinal momentum 
fraction u before the parton emits soft gluons. Of course, gp is not physical since we cannot 
separate a collinear parton from a cloud of soft gluons. Only after gp{u) is combined with 
the effect of the soft gluon emission, the parton distribution fp is physically meaningful. 
The relation between fp and gp is given by 

k)= JdujJ ^e'(-/2-?/"-P). J d7^e'^''-^>{n ■ P)5(r/) J du6{uj - 2un ■ P)g%{u). (38) 

Note that uj and yn ■ P are the label momenta, and 77 are the residual momenta. Therefore 
using the fact that 



and integrating the delta function with respect to uj yield 



(39) 



rp{y,K) = {n-P)S{K)gl,{y). 



(40) 
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In terms of the parton distribution function fp{y, k), T^^^ is written as 

t/°) = -C\Q) JdooJdn-kJ^J dr]e'^^/^+''''-''''-'^>Jp{n-k){n-P)S{T]) 
X J du5{uj — 2un ■ P)gp{u) 

= —C'^{Q) j dn - k j drj j duJp{n ■ k){n- P)S{r])S{un - P + n- q — n-k — i])gp{u) 
= -C^Q) Jdr]J duJpiun ■ P + n ■ q - 7]){n ■ P)S{r))g%{u) 

= -C^iQ) Jd7]J duJpiun -P + n-q- r])fp{u, r]). (41) 

In deriving this result, note that un-P is the label momentum of the parton, that is, n-p, and 
the exponent indicates the momentum conservation since a slight mismatch between 2un - P 
and the photon momentum n - q gives n - k + rj. The forward scattering amplitude is given 
by a double convolution of the jet function with the coUinear matrix element and the soft 
Wilson line. Because the jet function is affected by both the coUinear momentum and the 
soft momentum, it is impossible to write T^^^ as a single convolution with the conventional 
parton distribution function fp{y) without the effect of the soft gluon emission. However, 
as will be shown below, the moment of fp is given by the product of the moment of the soft 
Wilson line and Qp. 

Away from the endpoint region, we can neglect the soft momentum rj, and in this limit 
the soft Wilson line cancels to give fp{y, 0) = gp{y) is the conventional parton distribution 
function. And we recover the result away from the endpoint region 

Ti"^ = -JdyH{Q,y)fUy), (42) 

where H{Q, y) is the hard function, which can be split into C^(Q) and the jet function near 
the endpoint region. The main difference is that the effect of the soft gluon emission cancels 
away from the endpoint region, whereas incomplete cancellation occurs near the endpoint 
region. This incomplete cancellation results in the presence of the soft Wilson line, which 
represents the real soft gluon emission in the process. 

VI. MOMENT ANALYSIS 

In order to consider Ti^^ in moment space, let us introduce rj = {1 — v)n - p. The jet 
function Jp{n - k) has support only for a positive argument, and from Eq. (HTl) . n - k is given 
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in terms of the partonic variable n ■ p as 

n ■ k = un ■ P + n ■ q — rj = n ■ p — wn • p — (1 — v)n ■ p = {v — w)n ■ p. (43) 
Therefore v should he w < v < 1, and T]-°^ is written in terms of the partonic variables as 
Ti''\x,Q) = -C\Q) f'dy^^ dv{n-pfsi{l-v)n-p)jpi{v~w)n-p). (44) 

We take the discontinuity of Ti'\x, Q) to obtain the structure function Fi{x, Q). The hard 
coefficient C{Q) and gp{y) are real, therefore the imaginary part arises from the product of 
S (j^l — v)n ■ p^ Jp(j^v — w)n-pj. The procedure of taking the discontinuity can be performed 
either in SCETj or in SCETn- Since Eq. f l44l) is the result obtained in SCETn, we describe 
how the imaginary part can be taken in SCETn- The discontinuity of T^^^^ comes from the 
jet function Jp only since the soft Wilson line is real due to the fact that it is hermitian, 
hence its vacuum expectation value is real. The detailed discussion of taking the imaginary 
part in SCETi and SCETn, and the proof that the imaginary parts in both theories are the 
same are presented in Appendix B. 

In SCETn, we obtain the flavor nonsinglet structure function as 

F,{x,Q) = C\Q)Y.e] f'%Uy) [' dv{n-pfs({l-v)n-p)—lmJp(iv-w)n-p), (45) 

where is the electric charge of the parton q. For simplicity, we omit the summation over 
the quark flavors q from now on. Note that J pin ■ k) is actually the propagator of the n- 
collinear fermion, so it is of the form l/(^ ■ k + iO^) modulo logarithms oin-k with radiative 
corrections. Let us define the dimensionless function Jp near the endpoint as 

n-p — \m Jp{v — w) = -Jp{—\ (46) 

and let us also define the dimensionless soft Wilson lines as 

~S{v) = n ■ pS({l - v)n ■ p) = I(0|tr5t5„5(l -v + !^)siSn\0). (47) 
\ ^ JS ^ n ■ p ^ 

Then Fi(x, Q), with the explicit renormalization scales, can be written as 

Fi(x,Q) =C2(Q,/io) f'%%{y,f^) /'-5(t;,/i)Jp(-,/io,/i), (48) 

Jx y Jw V V 

where /xq ~ QV^ ~ x is the scale between SCETj and SCETn, and /i is the renormalization 
scale in SCETn with Q ^ /xq ~ QV^ — ^ ^ f^- According to the two-step matching, the 
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Wilson coefficient obtained from the matching at Q evolves down to the scale /iq. The jet 
function is computed from the matching between SCETj and SCETn at yUo, and it evolves 
to the scale /x. The soft Wilson line and the coUinear matrix element are evaluated at the 
final scale yU. This is the result in SCET in comparison to the result obtained in the full 



QCD factorization approach 23 1. 
If we write 



B{w)= f'-S{v)Jp{-), (49) 

Jw V V 

with w = x/y, the moment of Fi{x, Q) can be written as 

Jo Jx y ^y^ 

= C\Q) dxx""^ C dy dw6{x-wy)gUy)B{w) 
Jo Jo Jo 

= C\Q) r dyy^-'gUy) [' dww--^B{w) = g%^^ ■ B^. (50) 
Jo Jo 

The n-th moment of B{w) can be written as 

Bn = C dww''-^Biw)= t dww''-^ C—Siv)Jp(-) 

Jo Jo Jw V ^ 

dww^~^ dz du5{w — uz)S{z)Jp[ — ] 

In In ^ 1) ' 



Jo Jo ' ' ' V ■ 

1 ^ rl 



= f dzZ^-^S{z) I duu''-^Jp{u) = Sn ■ Jpn- (51) 

Jo Jo 
Finally, the moment of the structure function is given as 

Fl,niQ) = Aio)^P,n(/iO, Ai) ■ 9p,nil^) ' 5'„(/i). (52) 

The parton distribution function fp{y, rj) in Eq. ( l40l) can be written in terms of y and v as 
fUy, v) = {n- p)S({l -v)n.py-^ = S{vf-^, (53) 

^ ^ y y 

and if we take the double moment of fp{y, v), it becomes 

rp,m,n = J dvv^-' J dyy-'S{v)^ = S^g'p^^_,. (54) 
In terms of the moment fpmni moment of the structure function is given by 

The advantage of SCET in obtaining Eq. fl52p is that each component can be computed 
independently using perturbation theory, and we can clearly understand how these terms 
arise in SCET. 
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(b) ^'^^^ ^ igTa-^^[S{uj - n ■ {p + p') - n ■ q) - S{u! - n ■ {p + p) + n ■ q)] 

FIG. 3: Feynman rules for the collinear operator O^. (a) the tree-level operator (b) the operator 
with a collinear gluon with incoming momentum q^^. 

VII. RADIATIVE CORRECTIONS 

We can compute the radiative corrections for each term in the factorized expression for 
Fi^n in Eq. fl52l) . Let us begin with the radiative correction for the coUinear part (^p. There is 
a dehcate point in computing the radiative correction of the coUinear part. In any collinear 
loop integral, we integrate over all the possible loop momentum and the loop momentum 
can reach the region in which collinear particles become soft. Since the collinear and the 
soft fields are regarded as distinct in SCET, we have to remove the soft contribution from 
the collinear part to avoid double counting. For this purpose, the zero-bin subtraction 



method is suggested jl9|. Whenever there is a collinear loop diagram, the loop integration 
is performed by counting the loop momentum as collinear. Then the integrand is rewritten 
by counting the loop momentum as soft, and the integral should be subtracted to include 
only the collinear contribution. Otherwise, when the soft contribution is included, the soft 
contribution is counted twice. This had been missing in SCET and was first pointed out 



by Ref. [19|. Some previous calculations are not affected by the zero-bin subtraction, but 
conceptually the zero-bin subtraction is the correct step to avoid double counting. DIS is 
one of the examples in which the zero-bin subtraction should be performed carefully. 
Let us define the collinear operator ,the matrix element of which yields gpiy), as 

0^^ = UWnp{u;-V^)WlU (56) 

The FejTiman rules for including a single gluon are shown in Fig. [3l And the Feynman 
diagrams for the radiative corrections of at one loop are shown in Fig. |H The naive 
radiative corrections without the zero-bin subtraction, using the dimensional regularization 
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with D = 4 — 2e, are given as 
Ma = Mb = -2ig^CF- 



d^l n-{l+p) 



6{uj — uj') — 5{uj — uo' — 2n ■ 1) 



asCFit'2 r 



47r 2e 



5{uj -uj') + — 



UJ fe{uj)9{uj' -uj) e{-uj)d{uj - uj') 



{uj - CJ') + 



Mr 



_ 2^ t f d^'l {D~2)ll 
^2 J (27r)^/2[(/+p)2]2 

asCFiLl2{u' - uj) 
47r 2e 



A2 



5{uj -uj' -2n-l) 



uj)e{uj' - uj) + 9{-uj)d{uj - uj') 



(57) 



where uj' = n ■ {p + p'). Note that the terms proportional to 6{uj) {uj > 0) in Eq. (1571) 
contribute to the quark distribution function, while those with uj < contribute to the 
antiquark distribution function. Therefore the sum of all the corrections contributing to the 
quark distribution function is given by 



27r 2e 



25{uj - uj') + 1 + _ ^^^(^y 



(58) 



{u' - u) , 

while the contribution to the antiquark distribution function is obtained by replacing uj and 
uj' by —UJ and —uj' respectively in Eq. flSHl) . 

The zero-bin contribution in each diagram is obtained by the loop integral in Eq. fl57j) . 
where the coUinear loop momentum covers the soft region in which n ■ / ~ A and n ■ / ~ 

d^l 1 



M 



^ ^2 7 (27r)^/2(^.;+p2/^.p)^./ 



5(0; — uj') — b{uj — uj' — 2n ■ I) 



d^l {D - 2)1] 



-,5{uj — uj' — 2n ■ I). 



(59) 



2(n-p)2y (27r)^/2[rr-/+pVn-p]2 

Here M° is suppressed by A^/Q^ and it becomes zero when performing the loop integration. 
The total zero-bin contribution becomes 



2 TV 

1 / —p 



(^-±)(^-ln 



e /Cj' — UJ\~l~t , 



/in ■ p 



—)5{uj-uj') 

■ n' 



V /AT), ■ nl V 9 / 



(cj' - OJ) 



(60) 




FIG. 4: Radiative corrections for the quark distribution operator at one loop. 
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Since 



we can write 



/oo roo \ 
driri~^~'^9{ri) = / drir]'^'" = 
-oo Jo eiTV €n 



v-'-'Oiv) = (— - —Uv) 



+ 



CuV ^iR 



(61) 



(62) 



^euv em' ^+ 

where the subscript means the "+" distribution. Using this relation, the terms proportional 
to (5(a; — a;') cancel, and the divergent part of the zero-bin contribution is written as 

2 TT euv — ^^)+ 

As will be shown below, this is exactly the same as the soft contribution from the radiative 
corrections for S{r]), and it should be subtracted from Eq. (1581) . 

The relation between the bare operator O^"^ and the renormalized operator can be 



written as 



Of{uj)= I dco'Zico,co')OfW) 



(64) 



where the counterterm Z{uj^uj') is given by 



agCf 
2ne 



-6iuj-uj')+ . / V 
2 [uj' — uj)+ 



f\2 



The renormalization group equation for O':^ is given by 



duj'^{uj,u')Of{u'), 



-e{j - uj)e{uj)\. (65) 



(66) 



where the anomalous dimension 7(cc', uo') is given by 



UsCf 



TT 



-5{u^ - u') 



-1 + (iu/uj') 



l\2 



(67) 



{uj' - UJ)+ 

In order to express Eq. ( l66l) in terms of dimensionless variables, let us write u = 2Ey, 
uj' = 2Ez, where E is the energy of the quark and < y, z < 1. The renormalization group 
equation Eq. (!66|) is written as 



OsCp f dz 



TT 



y^ , -l + {y/zf 



Pqq{x 



as dx 

7T Jy X 

where Pqq{x) is the quark splitting function 

r3 



2^^ ,J ' {1-y/zU 
2Cf 



Of{z) 



(l-x)_ 



Pqqi^) — 



2'^'-"^"'(r^ 



(68) 



(69) 
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Note that, in Eq. (l68|l . there is an additional term —2Cf/{^ — x)+ due to the zero-bin 
subtraction. Therefore the matrix element g], of scales differently from the conventional 
parton distribution function fp{y) away from the endpoint region. However, when we include 
the effects of the soft Wilson line, we obtain the same result as the conventional approach. 
(See below.) The moment Qp^ satisfies the renormalization group equation 

d 



^^-^9p,n = -7c ■ 9p,n, 



(70) 



where the anomalous dimension 7c is given as 



7(7 = — — / dxx 

71 Jo 



n-1 



p ( \ _ '^^F 
'"^''^ (l-x)_ 



271 



(41niV-3) -41nA^ 



The last expression is the large n limit with = ne"^^. The first parenthesis comes from 
the splitting function Pqq{x) and the second parenthesis comes from the zero-bin subtraction 
-2Cp/{l-x)+. 

We now turn to the radiative correction for the soft Wilson line S{ri). The radiative 
correction for the soft Wilson line was computed in Ref. |16|], and we quote the result. The 
relation between the bare operator Spiv) the renormalized operator Sp{r]) is given by 



Sp{v) = Jdv'Z^,^{v,v')SR{v'), 



(72) 



where 



a^Cp 1 e{r] - 7]') 



71 e{T]-r]')+' 

The renormalization group equation for the dimensionless S{v) is obtained by putting rj 
(1 — v)n ■ p, 1]' = {1 — v')n ■ p, and it is given as 

2asCp /•! dv' S{v/v') 



(73) 



d ~, , 2asCp /•! dv' Siv') 

li—biv) = / — 

d/j, 7r Jv v' h— v/v' 



71 Jv V' (1 — f ') + 



(74) 



The n-th moment of the soft Wilson line S{v) satisfies the renormalization group equation 

(75) 



fJ'-T^n = -Is ■ Sn, 
dfl 



and the anomalous dimension 75 is given as 

2asCp 



Is 



71 



i/n-i-^41niV, 



(76) 



where Hn = Z]j=i l/ji aiid the large n limit is taken in the final expression. 
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Note that 75 is exactly the zero-bin contribution, as can be seen in Eq. (17T|) . Then the 
double moment fp^n,n+i = ^ndp^n satisfies the renormalization equation 

d 



" pq 



(77) 



where the anomalous dimension 

asCp 



7n = 7c + 7s 



271 



1 



n[n + 



2 1 



271 



(41nA^-3) 



(78) 



is the one obtained from the Altarelli-Parisi kernel. Therefore in moment space, the (double) 
moment of the parton distribution function even in the endpoint region satisfies the same 
renormalization group equation away from the endpoint region, or the result in the full 
theory. 

Finally, let us consider the radiative corrections to the jet function. To order as, the jet 
function Jp{n ■ k) is given as 



Jp{n ■ k) 



1 



n ■ k + ie 



1 + 



Att 



(21n^ 



—Pn ■ k — ie 



- 31n- 



—Pn ■ k — ie 



(79) 



where P is the label momentum P = n ■ p' = Q. Therefore Jq( (f — w)n ■ pj is written as 



Jgiiv - w)n-p 



1 



[v — w)n ■ p + ie 



A71 



(21n^ 



— (f — w)Qn ■ p — ie 



-3 In 



— (f — w)Qn ■ p — ie 

7 



7-^ 
3 



(80) 



The imaginary part of Jq{(v — w)n ■ pj is given by 

Am.jQy[v — w)n ■ Pj = d[v — w) IH (2m 



TT 



Arc 



- 3 In ^ + 7-71^ 



+ 



a.Cf 



An 



1 



41n^-3) + 



A ln(f — w) 



(81) 



The dimensionless jet function Jq{w/v), defined in Eq. (I46l) . with n ■ p = Q/w can be 
written as 



Jo I™ 



+ 



5(1-^ 

^ V 
OtsCp 



l + ^(21n^^-31n^ + 7 

Att V /,2 ,,2 



An 



TT 



An: (1 — w/v)- 



41n-^ -3 + 41n(l -Wt;) 

/i2 



(82) 
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where we neglect the In v term as f ^ 1 . The moment of Jq is given as 
f-i 



jQ,n 



duu"^ ^Jo(u) 





1 + 
1 + 



asCp 



in 



21n2^-31n^ + 7-7r2 + (41n^-3)iJ„_i-^ 



Q 



n-l 



4m 



k=l 



(21n2 



31n^ + 7- 



(83) 



where the last expression is obtained in the large n limit, which is consistent with the result 
in Ref. 3. 

We can present the moment of Fi{x,Q) to order a^. It is the product of the square of 
the hard coefficient C (Q) , twice the running of the hard coefficient from Q to Q / \fN using 
Eq. fl2T|) . the jet function at Q/^/N, and the running from Q/^/N to yU using Eqs. flTTl) . fl76l) : 



FiAQ) = C\Q,Q/^N)Jp,4Q/^N,fi)SMgj,^M 

l + ^f-16 + -) + ^f-21n^iV + 61nAr 



An 



3 ;+^(7c + 7.)ln^j 



l + ^(21n2iV + 31niV 
47r ^ 



y - 9) + (7c + 7s) In ^ 



'§n(/^)^|>,.(/x). (84) 



This result can be compared to the DIS structure function to one loop in Ref. 2J], where 
the moments of the nonsinglet structure function F2/{2x) = F\ are given as 



N 



l + ^52^^ + 7,ln^jA^(/i). 



(85) 



Here An{^) are the matrix elements of the twist-two operators renormalized at /x, and 7^ is 
equal to 7n = 7c + 75, given in Eq. (1781) . In the large N limit, is given by 

^2 



21n2A^ + 31nA^ 



9 



(86) 



Eqs. (|84l) and ( |85l) are the same, which means that the result for the moment of the structure 
function in the full theory away from the endpoint region can be extended to the endpoint 
region. However, our result is obtained near the endpoint region where the effect of the soft 
gluon emission from the soft Wilson line is present. The fact that the full-theory result can 
be extended to the endpoint region results from the relation fp^n,n+i = ^nQp^n in moment 
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space. Here the soft gluon emission plays an important role, and the effect shows up not 
only in DIS, but also in other high-energy processes such as Drell-Yan pro cesses and jet 
production in e~^e~ collisions. In the full theory it was considered in Ref. 251 ]. 

Eq. flMl) is the result at order a^. We can use the renormalization group equation to 
sum up all the large logarithms. The moments of the structure function in SCET to leading 
logarithmic accuracy is given by 



FiAQ) = C\Q)e 



-2h(Q,Q/ 



N] ~ 



Q \ -h(Q/VN,p) X 



where 



When we resum the large logarithms in Eq. (IHHi) . it is written as 

'a (^g/^)n2CF(3-87r/{/3oa,(Q)))//3o 



OlsiQ) 



Jp,n{Q/\fN)e-''' = Jp,n{Q/\fN[ 



as{Q/VN) 



(87) 



(88) 



(89) 



with Pq = 11 — 2n//3. Eq. is obtained by expanding Eq. (157|) to first order in a^. 



VIII. FACTORIZATION OF Fl(x,Q) 



The longitudinal structure function Fl{x,Q) is proportional to the imaginary part of 
T2(x, Q), as in Eq. (ITB]) . If we consider the tensor structure of the time-ordered products of 
the currents in SCET, T2(x, Q) is obtained by the products of the subleading currents Jj^""^ 
and J^"). Since the tree-level amplitudes vanish, we consider the operators with n-coUinear 
gluons from Eqs. ([23]) and 

1 



j{la) 



-n, 



1 



n ■ V 



(90) 



Here the delta functions are included for convenience, and Ba-, Bi, are the Wilson coefficients 
for the subleading current operators. As explained in Section II, the possible four types of the 
time-ordered products with Jj^""^ and Jj^^^ should contribute in the same way due to current 
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FIG. 5: The Feynman diagram describing the forward scattering amplitude T2{x,Q). 



conservation. Here we choose the two possible time-ordered products T[jj^^\z) JI)'^\0)] and 

verify that both contributions are the same by explicit calculation. 
And we show that the expression for the longitudinal structure function also factorizes. 
From the time-ordered product with J^'^\z) and JI)^"'\0), T2{x,Q) is written as 



X T 



X 



n-V 



X T 



X 



1 



(91) 



where the last expression is obtained by factorizing the usoft interactions after redefining 
the coUinear fields. The Feynman diagram for T2{x,Q) is schematically shown in Fig. [51 

Since there are no collinear particles in the direction in the final state, we can define 
the jet function Jp{uj, n ■ k) as 

J— [wlUpiW^ [5{uj - n ■ V)Wlin] [z ' ' ^ 

n ■ P ^ J L J 



(o|r 



n 



(0) 



|0) 



2] {2t^Y 



Jp{uj,n- k), 



(92) 



and T^"' can be written, with z = n ■ z/2, as 



fr(x,Q) 



dujdujB^{uj, uj) 



dn ■ kdz 
2^ 



(93) 



X T 
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/■ , ,_n2/ -^ fdn-kdz f, i((l-w)u>/2-n-k-v)z 

- J dujdujB^{uj,uj) J — J dr]e v ' ' V Jp{uj^ 



n ■ k) 



X UQ\ir[slSJ{n + n ■ id)SiSn] |0) ■T[^^WJ{u - n ■ V+)^W^^ , (94) 



N 

By putting 77 = (1 — v)n ■ p {w < v < 1), the spin-averaged matrix element between the 
proton state, Tg" is given as 

T^''(x,Q) = -^-P j dwdwBl{uj,w) j driS(ri)J^(u;,{l-w)^-r]) 

X C dy5{u-2yn-P)g%{y) 
Jo 

= - f'—gUy) fdvS{v) j duBl{Q,u)n-prp{u,{v-w)n-p), (95) 

We can also compute the contribution to T2(x, Q) using J^^^\ which comes from the part 
proportional to 11^111, in T^i,. It is written as 

(a;,g) = i j d^ze'^'^-P'+P^' j dudu'dcJBb{uj,uJ)Bb{uj' ,uJ) 

X T InWn^ [W^^W^] [5{u; - n ■ V)WlCn] {z) 

X %WJ{u;' - n ■ V')] \wl4)iwMcJ - V+)^WlUQ) 
By defining the jet function Jp(a;, n • /c) as 



(96) 



(o|r 



\wl^iw:\ - n ■ V)Wlin] {z) ■ [inWn] [wlUpiWr] (0) 



|0) 



.7^ r d^k 



_ .1^ [ 
= '2J 



e-'^-'J'piuj.n-k), 



(27r)^ 



after factorizing the usoft interactions is given as 
f^\x,Q) = - I du;dwB',{u;,uJ) | ^^I^^e'0'-"^''/'-"-')^4(u;, n • A;) 



(97) 



(98) 



X T 



/■ J J— n2/ — \ fdn-kdz f i( {l-w)uJ/2-n-k-'n) z ^, 

= -y dijjdujB^{uj ,W) j — ^ j dr]e \^ ' ' V J] 

X ho\tr{siSJ{r) + n-zd)SiS^)\0)-l^W^ 



p{oj, n ■ k) 



:-o[uj — n ■ Fj, 



and the spin-averaged matrix element between the proton state, T2^, is given as 
T^\x, Q) = -n-P j dujdu;Bl{uj,uJ) J dr]S{r])J^p{uj, (1 - w)w/2 - r]) 

X ^ / dy5{uj- 2yn ■ P)gl{y) 

U] Jo 

= - t—g%{y) tdv~S{v) I du;^^^^^jUu;,{v-w)n-p). 



(99) 
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a, a 




P j{lb) p' 



a, a 




= g'^T''T''-^{ S{uj + n-qi-n-p') 



+S{u! — n ■ p') 



n ■ p 



[2±= 7±/J 



n- qi \ n- q2 



71 -gi n-q2 n ■ {qi + q-j) \n ■ q2 n ■ qi , 

FIG. 6: Feynman rules of the current j^^^\uj) with one or two n-colhnear gluons. 

We can clearly see that both T^'* and factorize. By taking the imaginary part of T^'* or 
T2'', the longitudinal structure function Fl{x,Q) can be written as 



Fl{x,Q) = 8 f'%Uy) C dv~S{v) f duBl{Q,u)n-p 



— IraJpiuo, {v — w)n ■ 

71 ^ 



P 



FLix,Q) = 8 / ^gUy) I dvSiv) f dcu^^^^ 
Jx y Jx J Ti ■ p 



TX 

^aa Cphh 



ImJ^(u;, (f (100) 



where the first (second) expression is obtained using T^"" (T2 

The jet function J\{uj^n ■ k) at order can be computed from Eq. (!97|) . The Feynman 
rules for j*^^^^ with one or two coUinear gluons in the direction are shown in Fig. [61 And 
the matrix element Mi, from Fig. [5] with Jj^''\ after extracting the overall factor l/{n ■ p)^, 
is given as 

« / (2^p-'^"*(" -"■('+")) IITTTF''" P 



-g^c/-iD-2)J^6{u-n.il + p') 
.ifiasCp 



-9{uj)e{n ■ p - uj) 1 + ln 



P{l + p'f 
uj{n ■ p' — uj){—n ■ k — ie) 



n ■ p'fi"^ 



(101) 



2 271 n-p' 

where we collect the finite terms only in the last expression. Here n-p' is the total coUinear 



—— -S(UJ - n - P ) 7±a - i±= 

n ■ [p ~ q) \ n ■ qJ 



FIG. 7: Feynman rules of the current j^^°'\uj) with a single n-collmear gluon. 
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momentum of the intermediate states, and can be replaced hy n ■ q = Q. Therefore the jet 
function Jp{uj,n ■ k) is given by 

-rh / ,N «sCfi^^/ n^/^ n r, , uj(Q — Lj)(—n ■ k — ie)-\ 
JUuj, n-k) = -^-9{cu)9{Q - + In ^ ^-^ ^] . (102) 

Comparing the definitions of Jp and Jp in Eqs. fl92l) and fl971) . Jp(co', n ■ /c) is given by 

J^(cj,n- fc) = -^4(c<;,n- A;). (103) 

It is easy to see this relation by looking at the Feynman rules for which are presented 
in Fig. [3 Note that there is an additional factor of l/{n ■ p')^ in the definition of Jp, and 
the nonzero contribution comes from the part proportional to 7^0,, of which the radiative 
corrections are the same using either j^^""^ or 

Let us simplify the expression for Fl{x,Q) in Eq. fllOOl) . At order as, since the jet 
functions are already at order as, the Wilson coefficients C{Q) and Ba^bi'^) take their tree- 
level values, that is, 1. Introducing the dimensionless variable r = ujQ, the jet function Jp 
is written as 

Jl({y - w)n ■ p) = ^^^r9{r)9{l - r){l + ln[-^r(l - r)(-(l - -) - ie)]}. (104) 



Therefore the imaginary part is given by 

Jpfr, -) = --lmj'p( (v - w)n ■ p] = ^^r9(r)9(l - r). (105) 
V t) / TT ^ ^ 27r 

At order as, the longitudinal structure function Fl{x,Q) becomes 

F,{x,Q) = 8^ f'%Uy)- ['dvS{v) = 8^ t-gUy) f dvw~S{v), (106) 
Ait Jx y w Jw Arc Jx y Jw 

where the first (second) expression comes from Tg" {T2')- The two expressions differ by 
a factor of w'^, which gives a subleading correction since 1 — w ~ K/Q. Therefore the 
contributions from Tg'*" and Tg'' to Fl{x, Q) are the same near the endpoint region at leading 
order in A/Q. In fact, the result of the electromagnetic current conservation goes further 
than the fact that T2 can be obtained using either jj^""^ and J^^^-*- It should hold to all orders 
in as, and the scaling behavior of the two currents should also be the same. We present the 
anomalous dimensions for J^^''^ at one loop in Appendix C. 
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The moment of Fl{x, Q) using the first expression in Eq. fllOOp is written as 

FlAQ) = C dxx^'-^FL^x.Q) 
Jo 

= — - — — / dxx^~^ / dygp{y) / dw6{x — wy)— / dvS{v) 
n Jo Jo Jo w Jw 

^"-^ f"dyy-'g%{y) [' dww^' ^vSiv) 



IT Jo Jo Jw V 



— — — ^ dyy^ ^g'piy) [ dww^ ^ [ du ( dv5{w — vu)vS{v) 



TC Jo Jo Jo Jo 



TT n — 1 Jo 



dyy'^-'gUy) / dvv--'S{v) 



TT n ' TT n 



where the last expression is obtained in the large n limit. Compared to the moment -Fi,n, 
Fi^n is suppressed by n, which confirms that the longitudinal structure function Fl{x, Q) is 
suppressed by A/Q compared to Fi{x,Q). In order to obtain the exponentiated form, we 
should compute the radiative corrections to next-to-leading order accuracy. This has not 
been done here, but all the logarithmic terms such as W n/n can be resummed from the 
factorization property in SCET, as suggested in the approach using the full theory |26|]. 



IX. CONCLUSION 



DIS near the endpoint region can be described in SCET. The factorization of the structure 
functions is explicitly shown to order as, and the moments of the structure functions are 
expressed as a product of the Wilson coefficients, the moments of the jet functions, the 
collinear matrix elements, and the soft Wilson lines. The radiative corrections for each 
component can be separately computed using perturbation theory. The structure function 
Fi{x,Q) starts at leading order in A, while the longitudinal structure function Fl{x,Q) 
starts from order A and ag- Therefore the Callan-Gross sum rule holds at leading order in 
A, and the corrections can be systematically computed in SCET. 

High-energy processes, such as DIS, Drell-Yan processes, hadron collisions, e~^e~ jets, 
can be described by SCET. And all these processes possess common features though the 
detailed dynamics are different. First, the scattering cross sections are factorized, and the 
short- distance physics and the long-distance physics are separated. There is a universal soft 
Wilson line describing soft gluon emissions near the endpoint region with the appropriate 
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prescription for the soft Wilson lines depending on the external collinear particles. And there 



is a contribution from the parton distribution 



unctions if the initial particles are hadrons. 



Compared to the approach in full QCD 27|, there is advantage in employing SCET in 



DIS near the endpoint region. First, the factorization property becomes transparent since 
SCET is formulated from the beginning to decouple the collinear and the soft degrees of 
freedom. And the power counting in powers of A can be systematically performed. Secondly, 
once the factorized form is given, each factor has a different physical origin and its radiative 
corrections and evolutions can be computed in perturbation theory. The hard coefficient 
C{Q,fio) comes from the hard physics of order Q and can be computed in matching the 
full theory and SCETj. The jet function Jp arises from the hard-collinear physics of order 
\/QA, and can be computed by matching SCETj and SCETn. The soft Wilson lines and the 
collinear matrix elements are combined to give the parton distribution functions. Since the 
collinear particles and the soft particles are decoupled, the radiative corrections for the soft 
Wilson line are governed only by the soft interactions in SCETn, while those for the collinear 
matrix elements are governed only by the collinear interactions. To guarantee this, and to 
avoid double counting, the zero-bin subtraction should be performed. Each contribution is 
clearly separated and SCET specifies the prescription for computing radiative corrections. 

In this paper we have considered the fiavor nonsinglet structure functions. To be complete, 
the fiavor singlet structure functions should be included. In this case we have to consider 
the collinear operators with gluons, which contribute to the gluon distribution function in 
the proton. Though it will be more involved because of the operator mixing, the procedure 
for showing the factorization is straightforward. The complete treatment of DIS near the 
endpoint region including the fiavor singlet structure function will be considered elsewhere. 
It will be interesting to see if other various high-energy processes near the endpoint region 
can have similar features as those in DIS. 

Note added: The original preprint version of this paper did not include the zero-bin 
contribution [19|], and the parton distribution function in SCET was not properly defined 
excluding the soft part. In this paper, the zero-bin subtraction is correctly performed to 
solve the double counting problem, and the parton distribution function is carefully defined 
in the endpoint region. Because of these, the results of the calculations and the conclusion 
of the paper have changed. However, in the meantime, several authors jl?! . [isl criticized this 
paper based on the original preprint version. We would like to comment on the criticism 
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which appeared in the hterature. 



In Ref. [l7|], the authors argue that there should be no extra soft contributions outside 
the parton distribution function, which is correct based on the original preprint. However, 
in this paper, we include the soft contributions as part of the parton distribution function, 
which also affects the jet function in the endpoint region, contrary to the case away from the 
endpoint. They also claim that the n-coUinear gluon exchange shown in Fig. [8] (b), (c) are 
not kinematically allowed. But all the n-collinear contributions in Fig. [8] should be included 
because the interaction of the collinear gluon with the quark happens inside the proton, and 
the resulting momentum of the quark undergoing the hard collision should be regarded as 
n collinear, not the momentum of the quark before it interacts with a collinear gluon. The 
detail is explained in Appendix A. 

The main criticism of Ref. 18| is that the double counting problem is not performed 
properly, and that the parton distribution function is not correctly identified if only the 
matrix elements of the collinear operators are included. In the current paper, the double 
counting is treated using the zero-bin subtraction method. For the parton distribution 
function, we include the effect of the soft gluons in the parton distribution function. Besides 
these, we agree with their claim that the structure function depends only on Q^, and (1 — 
x)Q'^, not on (1 — which is illustrated clearly in Eq. (HHI) . where the structure function 
depends only on and /Xq = (1 — x)Q'^, and the dependence on /i cancels. 
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APPENDIX A: ZERO-BIN SUBTRACTION BEFORE THE USOFT FACTOR- 
IZATION 

For the purpose of illustrating the zero-bin subtraction method in SCETj, we consider 
the radiative corrections to the forward scattering amplitude before the usoft factorization. 
The starting point is that Tj^^ from Eq. ( l27l) is given by 
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(a) (6) (c) 

FIG. 8: Feynman diagrams of the contribution from collinear gluons to the forward scattering 
amphtude in SCETj. The mirror images of (a) and (b) are omitted. 

= zC\Q) J d^ze^('^+^-^')-^T[e^iy^7MW^ien(^)LW^n7.W^ien(0)], (Al) 

where the colhnear field is not redefined by Eq. fl25|) . and can interact with usoft gluons. 

We consider the radiative corrections with the n-collinear gluons, which contribute to 
the renormalization of the quark distribution functions. The Feynman diagrams are shown 
in Fig. [HI We consider the amplitudes only at the parton level,and the convolution with 
the parton distribution function is straightforward. The naive contributions proportional to 
—g^u from Fig. [8] (a), (b) with their mirror images and (c) are given as 

t/l 1 f d^l n-{l + p) 



Ma 

Mr 



2n-p'J {2t()^ P{l+pYn-V 
ii f d^l n- {I +p) 

2 " 



(271)^ P{l + pyn- {I + p')n-V' 



2^ ^ 



(A2) 



2 J {2n)D P[{1 + p)'^]'^n ■ {I + p') 
where p'^ = Pfi + is the n-collinear momentum, and we put p_L = for simplicity. 

Note that in Mb, the propagators are written in such a way that I + p' is collinear in 
the n direction. This should be included in the collinear contribution, contrary to the claim 




FIG. 9: When an n gluon is exchanged, there is a kinematical region where the initial state with 
/ + p is n colhnear, and the intermediate state with I + p' \s n colhnear, while the loop momentum 
I is n cohinear. 
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in Ref. |17|, in which the authors claim that Fig. [8] (b), and (c) should not be included 
since they are kinematically forbidden. However, by looking into the kinematics carefully, 
there are coUinear contributions from Fig. [H] (b), and (c). The point is that a quark inside 
the proton can interact with rl-collinear gluons before the hard collision. Therefore the 
rl-collinear gluon is regarded as part of the proton, and forms an n collinear jet. The 
parton distribution function describes the partons which undergo a hard collision after all 
the interactions with the collinear jet. 

The situation is schematically shown in Fig. [9], which is Fig. [8] (c). The collinear quark 
interacts with a collinear gluon before it collides with a hard photon. Therefore the longi- 
tudinal momentum of the collinear quark for the hard collision is n ■ (/ + p) , not n ■ p, where 
n ■ p is the longitudinal momentum fraction before it interacts with a collinear gluon. In 
order to see if the rl-collinear gluon is allowed by kinematics, let us introduce the partonic 
variable w', to avoid confusion, which is given by 

w 77— ^- (A3) 

n- [I +p) 

And n ■ {I + p) and n ■ {I + p') are given by 

1 -w' 

n-{l + p) = Q/w', n-{l + p') = n-{l+p + q) = {l-w')n-{l + p) = —Q. (A4) 

w' 

In the endpoint region where w' ^ 1, l + p can be n collinear, l + p' can be n collinear, while 

/ is n collinear. Therefore the contribution from rl-collinear gluons should be included. 

The reason why Eq. ( ]A2I) is naive is because the loop momentum can be soft, which 

should be avoided in the collinear sector. Therefore we subtract the contribution where the 

loop momentum becomes soft, and we call this the zero-bin contribution. It can be obtained 

from Eq. (]A2|) by power counting, where all the components of the loop momentum scales 

as A. The zero-bin amplitudes are given as 

— A - '^r ^ 1 f d^l n ■ p 

" ^ ^2^^ J {2TT)^P{n-pn-l + p^)n-V 

M, - -Atg Cf-J ^2n)^P{n-pn-l+p^)n-{l+p')n-V ^ ' 

where we put p"^ to regulate the infrared divergence. And the zero-bin contribution from 

Fig. M (c) is suppressed, and we neglect it here. The total zero-bin contribution is given as 

^ f d^l n ■ p r 1 1 



(2-k)^ P{n ■ pn ■ I + p'^)n ■ l^n ■ p' n ■ {I + p') 

A- ^ ^ ( J-1 ^lll. (Kf\\ 

^2n-p' J {27T)Dp{n-pn-l+p^)n-{l+p')' ^ ' 
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The usoft contributions from Fig. [TO] are given by 

M^s = WCft; 7 / , , 9N TT—^' (^7) 



2n ■ p' J {27i)^ P{n ■ pn ■ I + p^)n ■ {I + p') 
which is exactly equal to the zero-bin contribution. Therefore the correct computation 
including the zero-bin subtraction becomes 

M„ + M, + M, - (M° + M°) + = Ma + M, + M„ (A8) 

which states that the naive coUinear contribution without the zero-bin subtraction gives the 
correct result. It is also true in SCETn after the soft factorization. The zero-bin contribution 
to the collinear operator is the same as the radiative correction to the soft Wilson line. 
In calculating Ma + Mf,, note that we can write Mf, as 

iji r d^l n-{l+p) 



{2n)D P {I +pyn- {I +p')n- 1 



= -A.n^r f J!^\^^:l + P} n-jl + p) 1 

^ ^2n-p' J {2'K)DVp{l + pfn-l P{1 + pfn ■ [l + p'Y' ^ ' 

where the first term is equal to —Ma. Therefore Ma + M^ is given as 

2 „ ?4 1 f d^l n-{l+p) 



Evaluating the n ■ I integral by contours, doing the 1_l integral, and using the substitution 
n ■ I = —zn ■ p gives the infinite part 



TT e2n-pJo {1 — w + iO'^){l — z — w + iO+) 

'l + ti?ln 77^)) (All) 



TT e2n ■ p \ — w ^ iO+ ^ —w -|- 

where n- p' = {\ — w)n ■ p. Similarly, Mc is given as 



271 e2n ■ p Jo 1 — w — z + iO+ 

■ ^^-L (l + (l-w;)n — ). (A12) 

271 e2n-p^ ^ ' + iO+ ^ ^ ' 




FIG. 10: Feynman diagram of the contribution from usoft gluons to the forward scattering ampli- 
tude in SCETi, and the mirror image is omitted. 
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And the wave function renormalization for the external quarks is given as 

asCp 11^ 1 



w.f . 



47r en-p21 — w + i^'^ 



(A13) 



Noting that 



Im ■ 



w 



1 - w + iiO+ 
we obtain 
1 



ln(w;-l-iO+) = -vr- 



w 



Im 



-7r5(l-M;), (A14) 



-Im {Ma + Mb + + M^.f. = 77 

vr In en ■ p 2 



1 + 



The tree- level amplitude is given by 

-^tree 

the imaginary part of which is 

1. 



1 i 1 



n ■ p21 — w + iO+ 



-Im Mtree = „ 

TT n ■ p 2 



-o(l — If). 



Adding Eqs. (lAlSP and (]A17p . we have 

1 1/l 



n ■ p2 



w] 



2ti e^2'' ' ' (l-w), 
from which the anomalous dimension is given as 



7 



(-5 1 - w + — j, 

71 ^2 1 — w 



(A15) 



(A16) 



(A17) 



(A18) 



(A19) 



(l-w). 

which is exactly the Altarelli-Parisi kernel. Note that this is the result including the zero- 
bin subtraction, and it corresponds to the radiative corrections for the sum of the collinear 
matrix element and the soft part. 



APPENDIX B: DISCONTINUITY IN SCETj AND SCETn 

It is possible to take the imaginary part of the forward scattering amplitude to obtain 
the structure function in SCETi as well as in SCETn. If we only consider the collinear 
interactions with the intermediate state, the computation produces the jet function and 
the discontinuity due to the collinear interactions is the same both in SCETj and SCETn. 
Therefore the issue here is how to take the discontinuity related to the (u)soft interactions. 
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We consider the (u)soft interactions with the intermediate state in both effective theories 
and show that the discontinuity is the same. Since we are interested in computing the 
anomalous dimension of the soft Wilson line, we focus on the ultraviolet divergent part. 

Let us consider the usoft interactions in SCETj and take the discontinuity. The relevant 
Feynman diagrams are shown in Fig. [TT] (a), where the curly lines are soft gluons. Using the 
dimensional regularization, the Feynman diagrams in Fig. [11] (a) are given as 

c/^Z 1 



Ma = -Aig^CpOnJ 



{2n)D (12 + ^0+) {n-{l+p) + zO+) (n ■ {I + p') + iO+) {n ■ p' + zO+) ' 

(Bl) 

where On is the operator 

On = InWnlf-lMin. (B2) 
Evaluating the n ■ I integral by contours, doing the 1_l integral give the infinite part 

Ma = ^ , r dn . I , ^ -. (B3) 

vre n-p' + iO+ J-oo n ■ {I + p' + iO+) ^ ' 

By putting n ■ p' = n ■ {p + q) = {1 — w)n ■ p and n ■ I = —zn ■ p, we obtain 

a^Cp On 1 /-^ , 1 asCp On ln(w - 1 - iO+) 

Ma = / dz = ^ B4 

vre n ■ pi — w + ^0+ Jo 1 — w — z + zO+ vre n-p 1 — w + ^0+ 

where we neglect the \nw term as w — > 1. The discontinuity in SCETi from the usoft 

interactions is given by 

-ImMa = ^ . (B5) 

TT vre (1 — w)+ n ■ p 



This analysis is similar to the analysis in Ref. [IJ], in which a single-step matching was 
performed. 

In SCETn, the Feynman diagram is shown in Fig. [TT] (b). It gives 

J [1 — w)n ■ p — 7] + Tie r]^ 





(b) 



FIG. 11: Feynman diagrams for computing the radiative corrections of the (u)soft Wilson line in 
(a) SCETi and (b) SCETn. 
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where the first term in the denominator is the coefficient (jet function at tree level) with 
the energy transfer rj to the soft gluon. The remaining part is the result of the soft loop 



calculation 



16|. Taking the imaginary part of Mf,, we have 



1. OisCp 1 On , . 

-lmMb = r , (B7) 

TT vre (1 — w)-^- n ■ p 

which is the same result as Eq. ( ]B5[) obtained in SCETj. 
APPENDIX C: ANOMALOUS DIMENSION OF J^"^ 

We present the calculation of the anomalous dimension for J^^^^ at one loop, and explain 
why it is the same for J^^"-'. The current J^^^' from Eq. (|90l) is given as 

= -'^M / dcoB.ico) [^nWnSico - n ■ V^)] [w^UpiW^] = -n,j^''\ (CI) 

where Bh{uj) is the Wilson coefficient, which is 1 at tree level. The Feynman rules for j^^'^^ 
is given in Fig. [6l and the Feynman diagrams for the radiative corrections are given in 
Fig. [121 In Fig. [T21 diagrams (a) to (e) are the radiative corrections from the ra-coUinear 
loop diagrams. Diagram (f) is from the n-coUinear loop diagram, and diagrams (g), (h) 
are the contributions from the soft loops. We employ the background gauge field method 
for the triple-gluon vertex, and use the dimensional regularization with D = A — 2e. The 
external momenta p^, p'^ and are kept to give infrared cutoff, and the poles in 1/e are of 
the ultraviolet origin. 

Considering the flow of momenta, p' — q is the total outgoing collinear momentum in the 
direction, and n ■ {p' — q) = Q is the large scale in DIS. We use the variables 

n ■ {p' — q) = Q, n ■ p' = uj' = Qv, n-q = uj' — Q = {1 — v)Q, uj = uQ, (C2) 

and the allowed kinematic regions in DIS is Q > uj > 0. We extract the terms proportional 
to 7^^, and the divergent terms from each category (n collinear, n collinear, and soft) using 
the dimensionless variables are given as 

iMn = i{Ma + Mb + + Md + Me) (C3) 

,2 1 



11 AT „2 1 „/: 

6{u - v){2Cf(- + -) - -(in ^ + lnM(l - m)1 + — In ■ ^ 



N(A-u-v , fuil-u-v) , 

— f ^ ^(1 — u)v ' 
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9{u — v) u 6{y — u) 



) 



iMn 



{u — f )+ V {v — u)^ 

1 f / M(l-M-f) 

Ne^^{l-u){l-v) 
iMf = T^Cf5{u - t;) (4 - - In ^ + - 



x]e{l-u-v) + ^ — ^^^^^—^e{u + V-1)} 



i{Mg + Mh) = T^6{u - v) 



e /i^ 

-2C.N {-p'){-q') 1 . {-p'){-p") 

— + — In = ^ ~ ^ 1^ = — — 

e n ■ pn ■ qn^ JSe n-pn-p'ii'^ 



where T ^ is the common factor, given as 



a 



V. = -^gTaQ 



7±M 



(C4) 



47r " n ■ p 

We can compute the above matrix elements using the zero-bin subtraction. The infrared 
poles in 1/ eiR cancel when we add the soft contributions and the zero-bin subtractions and 
all the remaining poles turn into the ultraviolet poles. This procedure is similar to the pullup 
mechanism in NRQCD |28|. 



The relation between the bare operator j^^^^ and the renormalized operator j^^"' is 



dvZB{u,v)j^j^''\v), 



(C5) 




FIG. 12: Feynman diagrams for the radiative corrections of j^^^^(a;) in SCETj at one loop, p is 
in the direction, and p', q are in the direction (g incoming). Diagrams (a) to (e) include 
n-collinear loop, (f) includes the n-collinear loop, and (g), (h) are the soft corrections. 
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where the operators are dimensionless operators expressed in terms of m, v instead of uj and 
uj'. The counterterm Zb including the wave function renormahzation is given by 



Zb{u,v) 



1 + 



47r 



- + In 



^-^ln(l-M) 

47r iVe ^ ' 



5{u — v) 



+ 



Air e 



■ {u — v)- 



+ 1 



u 



/U 1 

+(-? — 

^7; 7; — 7 



-f (f — m)+ 



-)e{u-v) 



+ 



47r A^e 



u 



1 — u — v) 



+ 1 



u^6{v — u) 



(C6) 



■u — t;) + 



Note that the mixture of the uhraviolet and infrared divergences such as [ln(— g^//i^)]/e 
in Eq. (1030 cancels when all the contributions are summed. The renormahzation group 
equation for the current operator j*^^^^ is written as 



djjL J 
where the anomalous dimension '~fB{u,v) is given as 



(C7) 



(C8) 



rasCp /3 

V9 



In 



IT 
TT 2 



a.. 1 



ln(l — u) 5{u — v) 



1 u 



+ 1 - M - 



-)e{u-v) 



u 



TT 2A^ 



u — V 



{l-u){l-v) 



u 



+ 1 



uj6{y — u) 



n - t;) + ^ -e{u + V-1) 



Those terms in Eq. (ICSP proportional to 5{u — v) come from all the contributions, but the 
remaining terms proportional to the theta functions originate from the n-collinear radiative 
corrections. Compared to the renormahzation of the subleading heavy-collinear currents in 
Ref . 129! , the contributions of the n-collinear radiative corrections are the same because the 
contributing Feynman diagrams are the same. But the soft and the n contributions should be 
different due to the difference of the back-to-back coUinear current and the heavy-to-collinear 
current. Specifically the contributions not proportional to 5{u — v) in our computation and 
in Eq. (ICS!) are the same. For J^^"'\ the radiative corrections can be obtained in the same 
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way as in the case of J(i^), and it satisfies Eq. ffTU^ . 
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